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ABSTRACT 

In models for gamma ray burst afterglows, it is normally assumed that the exter- 
nal shock strongly amplifies the magnetic field and that this field maintains a steady 
value throughout the shocked region. We discuss the effects of modifying this (prob- 
ably simplistic) assumption by allowing for a short decay time. The observations are 
incompatible with a post-shock field that decays too rapidly. However if the field per- 
vades only a few percent of the total thickness of the shocked shell (and the electrons 
undergo only inverse Compton losses in the remainder) the model could be compatible 
with data. This would suggest a strong dependence of model parameters on the uncer- 
tainties of shock physics, therefore calling for independent external density estimates. 
We claim that afterglow emission should be instead seen as a laboratory where we 
can understand relativistic shock physics. The model we propose here together with 
afterglow observations in all wavebands could help to pin down the field structure. 

Key words: Gamma-rays: burst — radiation mechanisms: non-thermal — magnetic 
field 



1 INTRODUCTION 

It is widely accepted that GRBs afterglow emission is due to 
non-thermal radiation mechanisms produced in relativistic 
shocks, as the ejected matter from a highly energetic explo- 
sion expands into the surrounding medium. With current fa- 
cilities we can observe at photon energies up to 10 keV and 
the modeling of the observed lightcurves shows in general a 
predominance of synchrotron emission over inverse Compton 
scattering on the synchrotron photons. (Panaitescu & Ku- 
mar 2001 (thereafter PK01)). Only in one case, GRB000926, 
has a Compton component been inferred (Harrison et al. 
2001). This fact sets an upper limit on the density of the 
medium surrounding the explosion, whose properties can 
help us to unveil GRBs progenitors. In the standard frame- 
work the dynamics of the remnant expansion follows an adi- 
abatic law. This behavior is expected when the electrons 
cannot cool on the expansion timescale or (even when they 
can) the bulk of the energy remains in relativistic protons 
and magnetic fields. The post shock magnetic field, neces- 
sary for the synchrotron emission, is mainly generated dur- 
ing the shock itself, because the shock compression of the 
pre-existing field alone would lead to a negligible magnetic 
energy per particle (e.g. Gruzinov 2001). The magnetic field 
is normally assumed to be uniform and its lengthscale to 
be of the order of the remnant scale. The electron energy 
distribution behind the expanding shock front is given by 
a power law, down to a minimum injection energy. Model- 



ing eight afterglow lightcurves within this framework has re- 
sulted in a wide range of density from ~ 10 -3 to ~ 27 cm -3 , 
strongly indicating a low density environment (PK01). This 
is inconsistent with other independent observations placing 
GRBs birth in a star forming region: SN bumps in after- 
glow lightcurves (e.g. Bloom et al. 1999, Lazzati et al 2001, 
Bloom et al. 2002), metal lines (Antonelli et al 2000, Piro 
et al 2000, Reeves et al 2002), afterglow localizations close 
to star forming regions (Bloom, Kulkarni, Djorgovski 2002, 
Jaunsen et al 2002) and estimations of high density environ- 
ment from the temporal evolution of external column density 
(Lazzati & Perna 2002). Another remarkable inference from 
this model is that the fraction of internal energy given to 
electrons (e e ) is always much larger then the fraction given 
to magnetic field (et,). This implies that the magnetic field 
is relatively weak and the Compton emission dominates the 
overall cooling of the electrons. 

Recently many publications have investigated the some- 
what poorly understood generation of magnetic field in 
strong relativistic collisionless shocks (e.g. Kazimura et al 
1998, Gruzinov and Waxman 1998, Medvedev & Loeb 1999, 
Gruzinov 2001). Kazimura et al (1998) found generation of 
a small-scale quasi-static magnetic field; similar results were 
found by Gruzinov (2001) at small simulation times but his 
longer run shows that the field quickly decays by Landau 
damping; as a result synchrotron emission could be insignif- 
icant. Medvedev & Loeb (1999) found more reassuring re- 
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suits. Relativistic two-stream magnetic instability can, they 
claim, generate stable randomly oriented strong magnetic 
field in the plane of the collisionless shock front, so syn- 
chrotron emission is possible. These contradictory results 
show our ignorance of the possible acceleration process, and 
of how (or if) a stable magnetic field is created there. In the 
standard model this ignorance is "parameterized " through 
e e and es: the simplistic assumption is made that a per- 
sistent magnetic field makes these parameters uniform in 
space and time. An observation that maybe suggests a time 
dependence of e e and es is the independent determination 
of these parameters in GRB 970508 afterglow by Wijers and 
Galama (1999) at 12 days and Frail, Waxman and Kulkarni 
(2000) at ~ 1 year: while the former found e e = 0.12 and 
es = 0.089 the latter inferred e e ~ es = 0.5. 

These uncertainties, coupled with the apparent prob- 
lems with the standard model, motivate us to explore the 
consequences for afterglow emission of a different scenario, 
where the propagation of the magnetic field to large scale 
is disfavored. In this contest we particularly focus on the 
possibility that external density could be larger, keeping the 
total energy to a reasonable value, without loosing a gen- 
eral agreement with data. In fact, at this stage, we deal 
with sparse set of measurements, almost never simultaneous 
at all wavelengths. These data can be modeled in different 
frameworks, (e.g. GRB 970508 (Frail et al 2000, Chevalier 
& Li 2000), GRB 980519 (Frail et al 2000b, Chevalier & 
Li 1999), GRB 000301C (Berger et al 2000, Li & Cheva- 
lier 2001), GRB 991208 (Galama et al 2000, Li & Chevalier 
2001 and GRB 000418 (Berger et al 2001)). This analysis is 
meant to investigate how robust are the current estimation 
of shock and density parameters. 

In this paper we present this more general theory for 
the GRB afterglow emission (§2); our conclusion (§3) is that 
a qualitative comparison of this model with data suggests 
the possibility that a shorter magnetic lengthscale can be 
involved, affecting the estimates of parameters, but only a 
quantitative modeling of data can give the conclusive an- 
swer. 



2 SHORT MAGNETIC LENGTHSCALE 
THEORY 

The assumption of a post shock magnetic field persisting 
during all the expansion time has three main consequences 
for the synchrotron (S) emission and spectrum. First the 
emitting region linear dimension is R ~ c x t exp , where t exp 
is the expansion time in the lab frame; this enters in the 
calculation of the observed peak flux F p . Then the observed 
cooling frequency v c is computed using the Lorenz factor, 
7c of the electrons that cool radiatively on a timescale equal 
to the remnant age. Finally the lengthscale for synchrotron 
self-absorption has to be compared with the dimension of 
the fireball in order to compute the observed absorption fre- 
quency v a ■ The only break frequency which is not affected by 
any assumption about the extension of the magnetic field, is 
the observed peak frequency v m , corresponding to the min- 
imum random Lorentz factor 7 m : it depends crucially only 
on the fractions of internal energy given to electrons, e e , and 
given to the magnetic field ei,. Throughout the entire emis- 
sion volume electrons cool also via inverse Compton (IC) on 
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Figure 1. A schematic picture of the emitting region (shocked 
ISM). In region A, just behind the shock front, there are magnetic 
field and newly accelerated electrons that emit synchrotron and 
SSC for a time A/c before downstreaming into the (unmagnetized) 
region B where they Compton scatter the S radiation coming from 
zone A. 



the synchrotron photons (Self Synchrotron Compton: SSC). 
This cooling process dominates the electrons cooling if the 
Compton parameter Y = Lic/Ls (IC and S luminosity ra- 
tio) is greater than one. This happens if et <S e e and current 
estimates find in most objects eb/e e ~ 10 _1 — 10 -2 , (Wijers 
& Galama 1999, Granot et al 1999, PK01). 

If the magnetic field behind the shock persists for an av- 
erage timescale of tj, < t exp , its lengthscale is A ~ c x tt> < R. 
In principle A could depend on time, however we do not 
have any consolidated understanding of this process, there- 
fore for simplicity we will assume throughout this paper that 
X/R — 8 is constant during the afterglow phase. At each 
time the emitting region is then divided in two zones (see 
Fig [l]) . The first zone (region A) is the layer with dimension 
A just behind the shock front, where there are newly accel- 
erated electrons and magnetic field: here electrons emit via 
S and SSC; the second zone (region B) of linear dimension 
R — A, is filled by electrons that have just cooled by S and 
SSC emissions and that now Compton scatter the radiation 
coming from the first layer. We suppose that the newly ac- 
celerated electrons are advected with the flow; they produce 
synchrotron radiation only for a short time (with a cooling 
break therefore only at very high energies) in region A near 
the shock front and then they emit most of the IC radiation 
while downstreaming through region B. On the other hand 
if the electrons diffused freely, so that they could move back 
into the magnetized region after having lost most of their 
energy via IC, then the synchrotron spectrum would have a 
lower-energy break. This is probably less likely therefore we 
concentrate on the first case. The total observed spectrum 
then results from the sum of region A and region B spectra. 



2.1 region A: S and SSC emission 

The consequences of A < R on the S spectrum in region A 
can be derived easily from what we mentioned before. The 
observed peak frequency v~ m remains unchanged, 

V~ m = fro, (1) 
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and 

- 1/2 
7m OC V m = 7„ 



(2) 



In our generalized model with 5 < 1 we denote quantities 
by the symbol 'q', where 'q' is defined for 5 = 1, (for the 
explicit definitions of quantities in the standard model we 
refer to Panaitescu & Kumar 2000). 

The comoving peak intensity is I' oc n'B'X' , where n', 
B' and A' are respectively the external density, the magnetic 
field and the emitting region linear dimension in the comov- 
ing frame; therefore the observed peak flux F p is smaller 
than F p by a factor 5: 

F p ~r 3 (h'B'\')^ = SxF p , (3) 

where F is the bulk Lorentz factor. Since the magnetic field 
persists for a shorter time, a smaller fraction of electrons cool 
significantly, before the magnetic field disappears in that 
region. The Lorentz factor of the electrons, whose cooling 
time is equal to the typical timescale of the system, "J c is 
therefore higher: 

7c oc - — = 7 C x y S~ x , (4) 

t' b B'Hl + Y) 7 KJ 



where y — (1 -\-Y)/(l + Y). The observed cooling frequency 
v c oc B' 7c 2 r is then 



v c = v c x y 2 S 2 . 



(•>) 



The synchrotron self-absorption optical thickness can 
be approximated as 

f a (u) cx ( \ (ylv~ v )-^* = f p (v/v p )~^ , (6) 

for v < lXp, where p = c when the electrons are radiative 
(7c < 7m) and p = m when they are adiabatic > 7m). 
We will refer to the first case as fast-cooling regime and 
to the second case as slow-cooling regime. The absorption 
frequency corresponds to T a (v a ) = 1 and in slow cooling 
regime is given by 

v a = v a X <5 3/5 , (7) 

where we use Eq. |l] and Eq. ^ in fast cooling regime Eqs |E| 
| and [jj (y = 1) give 

v a = v a x S s/5 . (8) 

Eqs |^, and ^ are calculated in the case v a < v p (see 
Granot & Sari 2001 for the different possible S spectra). 
Nevertheless v a > v p for high densities 



h > 7 x W Z S 



2 r-12/7 F -12/10 



(9) 



in fast cooling and in slow cooling 

n > 8.2 x \Q 5 5- l E^ 2 e\% e%_ 2 (1 + z) t 5 J 2 y (10) 

where E is the isotropic equivalent energy and td ay is the 
observed time in days. The inferred external density is 
reasonable for 5 = 0.1 (n > 3.6 x 10 4 cm -3 and h > 
8.5 x 10 6 t^ 2 ( cm" 3 respectively) and for 8 = 0.01 (n > 

2 x 10 6 cm" 3 and n > 8.5 x 10 7 (^cm" 3 J but for lower 
S the density lower limits expressed in Eq g and Eq [l^ are 
perhaps too high for the afterglow environment. For sim- 
plicity we focus thereafter on i> a < v c but the reader should 



keep in mind that these cases hold as long as Eq ^| and Eq |l(j 
are not satisfied. 

The magnetic energy -|^- oc nT 2 and the Compton pa- 
rameter are higher at the beginning of the afterglow evo- 
lution and so the electrons are more likely then to be ra- 
diative. Since the minimum injected y m oc F decreases with 
time while 7 C increases as the efficiency of radiative loss de- 
creases, the electrons undergo a transition to the adiabatic 
regime. The transition time between the fast cooling and the 
slow cooling regime is 



f fa = Tfs x 8 2 
for 



(11) 



ft < 40 5~ 2 E^ 1 (Y F + I)" 2 e^"-! (1 + z)- 1 t~ d l y (12) 

the slow-cooling regime dominates the fireball evolution dur- 
ing times when the afterglow is observed ( ^ 1 day). 

The accelerated electrons that upscatter synchrotron 
photons give rise to a SSC component in the spectrum (for 
the SSC in the standard model see Sari & Esin 2000). Dur- 
ing the afterglow phase the medium has a Thomson optical 
depth 



3 



(13) 



of the order of f ~ 10 -5 x 5 for h = 10 2 cm -3 (a T is the 
Thomson cross section), therefore only a minor fraction of 
the emission is upscattered, and the synchrotron flux is not 
significantly altered. Nevertheless if Y > 1 the electron cool- 
ing is Compton dominated and SSC upscatters low energy 
photons, that might be detected as a Compton component 
in x-ray and, in the future, in higher energy band spectra 
(Sari, Narayan & Piran 1996). The parameter Y can be de- 
fined as the mean number of scatterings times the average 
fractional energy change per scattering, gained by a photon 
traversing a finite medium (Rybicki & Lihgtman 1979). For 
a Thomson thin, relativistic medium Y is 



(14) 



where N^) is the normalized ~y distribution (see e.g. PK00 
for explicit formulae for u a < v p ). For Y > 1 one may need to 
consider more than one scattering per photon; however the 
energy of a photon in the rest frame of the second scattering 
electron is a factor 7 higher than the original energy, gen- 
erally exceeding the Thompson regime limit (511 keV); the 
scattering cross section is then substantially reduced (Klein- 
Nishina (KN) regime) and further scatterings inhibited. As 
a consequence we will consider only single scattering of S 
photons. The KN effect limits also the electrons that mostly 
contribute to the SSC luminosity. They have a Lorentz fac- 
tor 



7 < Ikn — 511 T/hv Ep , 



(15) 



where ue p is the S energy peak frequency (y c in slow cool- 
ing and Vm in fast cooling). In zone A the S spectrum can 
have a very high v c (see Eq|B|); consequently in slow cooling 
7c could be greater than ^kn and the Klein-Nishina effect 
substantially decreases Y (we will refer to this situation as 
"zone A in KN regime"). Moreover if Y ~ 1 (which is the 
case in most simulations) we can omit the contribution from 
7 > 7xat. Eq. ^ for 2 < p < 3 gives then 
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Y s (l + Y s ) {3 - p) =Ys(l + Y s ) {3 - p) x S (p ~ 2 \ (16) 
for % < 7k jv and 

ys(i + ys) (2(p ~ 3)) = is(i + y s ) (2(p_3)) xs (7 - 3p \ (n) 

for 7c < 7k jv. If the electrons are in the fast cooling regime 
it is very unlikely that 7™ > 7k at and Eq. |l4| gives 



Yf = g -r 7m 7 c = Y F ; 



(18) 



note that Yf does not depend on 8, because most of the 
electrons cool before the field disappears. Given a S seed 
spectrum Fs(v s ), the resulting single scattering IC spectrum 
is (Rybicki & Lightman, 1979) 



fJ C =A I F S 



4-y 2 X 



f(X) dX 



N(y) dry, (19) 



where A = 3arR8 (in this model), X = u/i^Us and 
f(X) = 2XlnX + X - 2X + 1 includes the exact cross 
section angular dependence in the limit 7>1 (Blumenthal 
and Gold, 1970). For the cases we are treating it is possible 
to use for the first integral the analytic expressions given by 
Sari & Esin (2000) in their Appendix A. The IC spectral 
peak flux and break frequencies for 8 < 1 are related to the 
IC spectrum for 8 = 1 in the following way 



F P W ~ tF v = Fl c x S 2 
using Eq |^ and Eq ^| 

— ^7m Vjxi 

using Eq. § § 



-IC ^ 9 - 2 ~ _ IC 
V m — £]m Vm — Z^ m , 



-IC o ~ 2 - IC , 

J/ c ~ 27c V c = U c X 



(20) 
(21) 
(22) 
(23) 



using Eq. g |; 

-IC „ o - 2 - /C „ c3/5 

in the slow cooling regime using Eq. ^[ |^ and 

-IC o - 2 - TC7 „ r- 2/5 / <n 

i/ a ~ 27c v a = u a x 8 1 (24) 
in the fast cooling regime, using Eq. ^. 

2.2 Region B: IC emission 

Zone B extends over a distance (R — A); since we are mainly 
interested in 8 < 10 , we approximate (R — A) ~ R 
throughout this paper. This region is populated by electrons 
that have radiatively cooled via S and SSC for a time tb. If 
they were originally injected at the front shock with a power- 
low energy distribution ^(7) oc 7 _p then the steady state 
distribution in region B depends on whether the electrons 
were adiabatic or radiative in region A. The 7 c b for electrons 
that cool in a dynamical timescale t exp is in fact 



7cs = 



i„ U rad + 
-<T T —h^- t dy 



(25) 



where U ra d is the S radiation energy density and the factor 
1 /2 takes into account that only half of the S photons travel 
from region A towards region B in the electrons rest frame. 
The calculation of U ra d in the two different regimes leads to 



7cs oc 



(n 2 etrf 
(n 2 etTym 



' 1 7c 3 " p R 2 sy 

%R 2 S)- 1 - 



7m < 7c 
Jm > 7c 



(26) 



If electrons in region A are in the slow cooling regime, the 
new steady state configuration is 



y P 7m < 7 < JcB 

N(j) oc t j-(P+ 1 1 7cS < 7 < 7c 



"(p+2) 



7 > 7c 



for slow cooling (7™ < 7 c s) or 



JV( 7 ) oc { 

-(P+2) 



7 2 7cs < 7 < 7" 



7m < 7 < 7c 
7 > 7c 



(27) 



(28) 



for fast cooling (7m > 7cs). 7cS > 7c is very unlikely, even 
if the electrons in region A have two sources of energy loss 
(S and SSC), because tb < t exp and SSC can be inhibited 
further for KN effect. If electrons in region A undergo fast 
cooling, the most common case is 7cs < 7c and the new 
steady state configuration is 



7 2 7cs < 7 < 7c 

iV(7) oc \ 7~ 3 7c < 7 < 7m 

7" (p+2) 7 > 7m 



(29) 



The "external Compton" spectrum in zone B is calcu- 
lated by means of Eq|l| with A — \(Jt R and the Compton 
parameter Yb using Eq [m] and Eq |2(] with tb = \- 
Even if zone B has a higher number of IC emitting elec- 
trons , the Compton parameter can be lower than in region 
A: a higher fraction of the electrons in region B cool effi- 
ciently and they are in average less energetic then in region 
A, therefore increasing the density parameter Yb decreases 
while Y increases because in zone A 7c moves towards low 
values and the KN effect is less and less important (until 
the electrons become radiative and Yf does not depend on 
density). In this situation the SSC radiation can actually 
account for most of the energy that goes into IC luminosity, 
(this is the case in Figs ^ and ^ right panels) ; anyway the 
SSC usually dominates the spectrum at high energy because 
it has a higher energy break. The peak flux is 



t b F p , 



(30) 



which is a factor 8~ l /2 higher than F p in region A; while 
the spectral breaks are related to the S seed spectrum in 
slow cooling as follows 



~ic 0-2 - 
v a -27 m ^ a , 

-IC ^ cy~2 - 
V m — ^7m V m , 



~IC ^, o-2 - 
V c ^ 27c B V c , 

if region B is in the slow cooling regime as well, and 

- ic 0-2 - 

Va ~27c B ^a, 
-IC r,~2 - 

v c - 27 m v c , 
v^ ^ 27c B i> m , 



(31) 
(32) 
(33) 

(34) 
(35) 
(36) 



if B is in the fast cooling regime. If S seed spectrum is in 
the fast cooling regime and 7 c s < 7c < 7m 



- ic ^ r,~2 - 

Va ^ 27c S V a , 

-IC r,~2 - 
V c — 27cB V Cl 

-IC n-2 - 
V m — 27 c V m . 



(37) 
(38) 
(39) 
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3 DISCUSSION AND CONCLUSIONS 

In Fig |^ and Fig. ^| left panels, we show as examples how 
the spectrum (respectively for fast and slow cooling) mod- 
ifies if a short magnetic lengthscale is taken into account. 
Note that S and IC peak fluxes are lower (see also Eq. || 
Eq. Eq. p(i| ), a lower energy break corresponding to the 
absorption frequency is present (see Eqs. |^ and ^ and the 
cooling frequency is at higher energies (see Eq. ^|). The SSC 
component in the slow cooling spectra of region A (Fig ^ 
is depleted because the KN effect is important. In the right 
panels we raise the density until F p cx h 1 ^ 2 5 = F p cx n 1//2 , 
thus n = w. We obtain v c = v c and the 8 < 1 spectrum 
fails to match the standard spectrum only at radio wave- 
lengths. In this band such high densities, in fact, overcom- 
pensate the shorter magnetic lengthscale and v a becomes 
greater then v a - Another important feature is that the total 
(region A plus region B) Y parameter is never much higher 
then Y, even when Y corresponds to a density 5 2 lower; in 
fact when F p — F p , the main contribution to IC luminosity 
comes from region A and Eq ^ and Eq [Hi] hold: only Ys 
depends on density and 5 but the dependence are very weak 
(Ys (1 + Ys) (3_p) oc n (p " 2)/2 x <5 (p ~ 2) ). These illustrative ex- 
amples of how the model works should now be followed by 
a modeling of real data. In fact the general agreement be- 
tween the two spectral shapes suggests that data could also 
be well reproduced in this framework, with possibly differ- 
ent values for the break frequencies and peak flux respect 
to the standard model. As regard the corresponding esti- 
mated density compared to the standard model, there are 
two competitive facts: the peak flux is brought to observed 
ranges by an higher density but the absorption frequency 
cannot be too high (see fig. 2 and 3) (note that unlike the 
standard model a high density, £ 10 2 cm -3 does not pro- 
duce an extreme inverse Compton component, which is not 
observed); Therefore, even if there is the possibility that a 
higher density is required, the suitable value can be derived 
within the errors only by a broadband modeling of data, 
with a consistent variations of all the parameters. 

Similarly to the standard model there are univocal cor- 
relations between the estimated spectral parameters and the 
derived fireball energy, external density and shock parame- 



ters. They are 

E cx i? p 3/ %- m - 5/ V c 1/4 /7 a - 5/6 ,5- 1/2 r a5 (40) 

n cx F p " 3/ Vn 25/12 ^ 3/4 ^ 25/6 5 1/2 T- 7/2 (41) 

i e cx F p - 1/2 ^ 12 ^ 4 u a 5/6 8^ 2 T-^ 2 (42) 

e~ B oc F p XI \~ m -^v-^v a -^H- z l % T h l\ (43) 
in slow cooling regime while for a radiative regime 

E cx F p 3/2 v c - 1/6 v a - 5/6 S- 1/2 T-°- 5 (44) 

n cx F p - 3/2 u c 17/6 v 25/6 8 1/2 T 7/2 (45) 

U CX F p - 1/2 ^ 1 /V c 2 /3 l7a S/6 5 l/2 r 3/2 (46) 

eh cx F p 1/2 v c - 5/4 V ' a - 5/2 5- 3/2 T- 5/2 . (47) 



Some constraints must be added to this set of equations. 
Because e e + es + e p = 1, where e p is the fractional energy 
that goes to protons, we most likely expect e e + en < 0.5. 



The kinetic energy in the afterglow should not be much 
lower then the energy released in the 7-ray phase, other- 
wise very large efficiencies of radiation, incompatible with 
internal shock models, would be required (Lazzati, Ghis- 
ellini & Celotti 1999). Finally the apparent source size, the 
rate of expansion and the energy in each afterglow should 
be consistent with observations in the radio wavebands of 
interstellar scintillation and its quenching (Waxman, Kulka- 
rni & Frail 1998). In principle these requirements together 
with observations of the spectral parameters at a particular 
time would set limits on 8 and therefore on the structure 
of the magnetic field generated in shocks. Unfortunately the 
current data do not allow a model independent estimation 
of spectral parameters. Actually the spectral shape, and in 
particular the break frequencies and the peak flux, are not 
unequivocally constrained by data in even the best studied 
afterglows (e.g. GRB 970508, GRB 000418). These spectral 
features are derived by modeling the sparse measurement 
collected in various wavebands at different times: very often 
different versions of the standard model can fit the data, 
consequently finding different spectral parameters (see for 
example Berger et al 2000, for GRB 000418). Moreover these 
model dependent quantities (above all i> a ) enter with higher 
dependences than S in the estimation of the fireball param- 
eters (Eq. po[]4^ ); therefore the break frequencies and peak 
flux, derived under the assumption that 5 = 1, cannot obvi- 
ously be used to figure out the value of S. These equations 
and the physical limits on the fireball parameters should be 
used instead to constrain the overall modeling of data. The 
results would actually lead to very interesting conclusions. 

If the modeling succeeded in reproducing the data with 
5 < 1, a different value for the external density, energy and 
shock parameters would possibly be derived. If it indicates a 
denser environment, this could be more compatible with the 
strong observational evidence that would place GRBs explo- 
sions in a star forming region. However a very short scale 
magnetic field (say 5 ^ 10~ 4 ) would be probably ruled out, 
since extremely high densities would be required in order to 
get fluxes in the observed ranges and the radio flux [y < v a ) 
would probably be too heavily absorbed: F v oc n -1 / 2 (in adi- 
abatic) and F v oc <5 _1 n~ 5 / 3 (in radiative regime). An other 
implication of the possibility to fit the very same data with 
models based on different shock physics assumptions is that 
parameter estimations greatly suffer from shock theory un- 
certainties. In this work we have only explored the effect of 
a short-scale magnetic field, but other reasonable variations 
of the standard model (such as time dependent eB,e e and 
p) may affect the parameter estimates as well. This would 
call for alternative model independent external density and 
energy estimates in order to constrain, together with the 
afterglow modeling, the shock physics. 

On the other hand if this model can reproduce data 
only with 5 = 1, it could be either evidence in favor of 
an extremely stable field or the indication that some other 
modifications of the theory must be taken into account. For 
example, as already mentioned, 5 can vary with time: for in- 
stance if the Rayleigh- Taylor instability occurs at the con- 
tact discontinuity, the lengthscale would grow as the fire- 
ball evolves and consequently late multibands observations 
would be more consistent with 5 = 1, while early measure- 
ments would require 5 <C 1. Observational tests of this pos- 
sibility can be performed using the theory developed in this 
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work; in fact we can look for a coherent changing in time 
of 6, modeling instantaneous spectra at different times, with 
all the remaining parameters kept constant. Of, course, with 
a very good dataset, this approach can be extended to all 
the other shock parameters, including es,e e and p. Such a 
procedure would yield unbiased estimates of shock param- 
eters variability, and the GRB afterglow would become a 
laboratory where shock physic can be tested. The best case 
up to now for constraining 8 and its temporal behavior is 
probably GRB 000926, but in a year SWIFT, REM and 
and possibly the VLA will provide us with more simultane- 
ous data (also at early times) in X-ray, optical, infrared and 
longer wavelengths and therefore more constraining fits will 
be performed. This work actually highlights the importance 
of observing afterglows at all wavelengths, from radio to hard 
x-ray, but emphasizes as well that current estimations from 
this data of the external density surrounding 7-ray bursts 
may depend on uncertain details of shock physics. 
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Figure 2. An example of how the spectrum modifies if S < 1. Total (solid line) and SSC+IC (dot line) spectra for <5 = 0.01 and total 
(3 dots-dash line) and SSC (dot-dash line) spectra for S = 1 at 1 day. In the left panel the two total spectra have all parameters, except 
for 8, the same: n=1.35 cm -3 , E=3 X 10 52 erg, e e = 6x 10 -2 , eg = 4 X 10~ 3 , p=2.2 at a redshift z=l. These parameters are typical 
of broadband afterglow lightcurves modeled in the standard framework (PK01). Both spectra are in the slow cooling regime (also zone 
B with -y c g < 7c). The SSC in zone A is highly inhibited by the KN effect; therefore most of the IC component is due to "external 
Compton" in zone B. The resulting Compton parameters are Y = 0.40 and Y = 1.05. In the right panel we raise the density in the 
S = 0.01 spectrum until F p = Fp, h = 1.35 X 10 4 cm -3 . Zone B is now in the fast cooling regime therefore ^ c g moves towards low values 
but since in zone A the electrons are almost all in Thompson regime the total Y increases, Y ~ 1.10: region A now dominates the IC 
emission. It is worth noticing that the two total spectra are distinguishable only in the radio waveband where the S = 0.01 flux is more 
absorbed. 
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Figure 3. Same as Fig. g with parameters n = 27 cm -3 , E = 1.8 x 10 53 erg,e e = 0.3 and eg = 8 X 10~ 3 and p=2.43 and S = 0.1. This 
set of parameters results from the best fit of 000926 (Harrison et al. 2001). The redshift of this burst is z = 2.0639. In the left panel both 
spectra are in fast cooling (also zone B with 7 c s < 7 C ). The Compton parameters are Y = 2.59 and Yp = 2.34. In the right panel we 
rise the value of density in the 5 = 0.1 spectrum until F p = F p , n = 2.7 X 10 3 cm -3 . The spectrum is obviously still in fast cooling with 
a lower f c g', therefore the total Y is lower , Y ~ 2.37 and the main contribution to the IC luminosity comes from region A. Again the 
two total spectra are distinguishable only in the radio waveband. 



